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We  investigate  numerically  the  electric  and  the  magnetic  response  of  the  split-ring  resonators  (SRRs)  for 
different  electromagnetic  (EM)  field  polarization  and  propagation  directions.  We  have  studied  the  current 
density  J  of  circular  and  rectangular  SRRs  at  THz  frequencies.  At  low  frequencies,  J  is  confined  to  the 
edges  of  the  SRRs,  while  at  high  frequencies  J  is  all  over  the  width  of  the  wire.  The  diamagnetic  response 
of  the  SRRs  is  also  examined.  Finally,  the  role  of  losses  and  the  magnetic  resonance  frequency  are  studied 
as  the  size  of  the  SRRs  becomes  less  than  100  nm. 
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1  Introduction 

The  idea  of  left-handed  (LH)  materials,  i.e.  materials  with  both  negative  electrical  permittivity  (s)  and 
magnetic  permeability  (//),  where  the  electric  field  (E),  magnetic  field  (B),  and  wave  vector  (k)  form  a 
left-handed  coordinate  system,  was  developed  by  Veselago  [1]  decades  ago.  However,  it  was  only  re¬ 
cently  that  such  materials  were  investigated  experimentally  [2-4],  despite  their  rich  physics  and  the 
large  number  of  associated  potential  applications,  e.g.  in  the  formation  of  a  perfect  lens  [5].  The  first 
experimental  realization  of  LH-materials  was  achieved  by  separately  constructing  ^  <  0  and  ju  <0  com¬ 
ponents,  and  then  by  combining  them  together  forming  a  left-handed  material  (LHM).  Although  it  has 
been  well  known  how  to  obtain  a  ^  <  0  material  easily  (e.g.  using  wire  arrays  [6]),  the  realization  of 
ju<0  response  (especially  at  high  frequencies)  has  been  a  challenge,  due  to  the  absence  of  naturally 
occurring  magnetic  materials  with  negative  jU .  The  possibility  of  the  realization  of  a  ju<0  material  was 
predicted  in  1999  by  Pendry  et  al.,  who  suggested  a  design  made  of  concentric  metallic  ring  with  gaps, 
called  split  ring  resonators  (SRRs)  which  exhibits  a  //  <  0  regime  in  the  vicinity  of  the  magnetic  reso¬ 
nance  frequency  co^  of  the  SRR  structure  [7]. 

In  the  past  few  years  there  has  been  ample  proof  for  the  existence  of  Negative  Index  Materials  (NIMs) 
in  the  GHz  frequency  range.  A  lot  of  groups  were  able  to  fabricate  [2,  8-10]  NIMs  with  an  index  of 
refraction  n  =  -\  with  losses  of  less  than  1  dB/cm.  Recently  different  groups  observed  indirectly  [11-14] 
negative  ju  at  the  THz  region.  In  most  of  the  THz  experiments  [11,  13,  14]  only  one  layer  of  SRRs  were 
fabricated  on  a  substrate  and  the  transmission,  T,  was  measured  only  for  propagation  perpendicular  to 


*  Corresponding  author:  e-mail:  soukoulis@anneslab.gov 

I^J^WILEY 

V  interScience® 

>t7-'  DI  =  COVf«  SO^.EIM,N5  SBE*r 


©  2007  WILEY- VCH  Verlag  GmbH  &  Co.  KGaA,  Weinheim 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

2007  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Magnetic  response  of  split  ring  resonators  at  terahertz  frequencies 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Ames  Laboratory  and  Department  of  Physics  and  Astronomy, Iowa  State 
University,  Ames,  Iowa  500II 

8.  PEREORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIEICATION  OE:  17.  LIMITATION  OE 

APl9!TR  apt 

18.  NUMBER  19a.  NAME  OE 

rtln  DAPIhQ  DTnQPrtMQTTiT  U  PlnDQ/TM 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  jj-jj 

unclassified  unclassified  unclassified 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


C.  M.  Soukoulis  et  al.:  Magnetic  response  of  split  ring  resonators  at  THz  frequencies 


1182 


the  plane  of  the  SRRs,  exploiting  the  coupling  of  the  electric  field  to  the  magnetic  resonance  of  the  SRR 
via  asymmetry  [15].  This  way  it  is  not  possible  to  drive  the  magnetic  permeability  negative.  Also,  no 
negative  n  with  small  imaginary  part  has  been  observed  yet  at  the  THz  region  [16,  17].  One  reason  is 
that  is  very  difficult  to  measure  with  the  existing  topology  of  SRRs  and  continuous  wires  both  the  trans¬ 
mission,  r,  and  reflection,  R,  along  the  direction  parallel  to  the  plane  of  the  SRRs.  So  there  is  a  need  for 
alternative,  improved  and  simplified  designs  that  can  be  easily  fabricated  and  experimentally  character¬ 
ized.  This  was  recently  achieved  in  the  GHz  region  [18,  19]  and  in  the  THz  region  [16,  17]  by  the  use  of 
finite  length  wires  and  the  fish-net  topology. 

In  this  manuscript,  we  will  present  results  of  the  magnetic  response  of  SRRs  as  the  size  of  SRRs 
reduce  in  size.  As  shown  Zhou  et.  al.  [20],  scaling  {co^  vs  ML)  is  obeyed  for  T>500nm,  while  for 
smaller  Z,  the  magnetic  resonance  frequency  saturates.  We  have  systematically  studied  the  current  den¬ 
sity  as  well  as  the  magnetic  response  of  the  SRRs.  We  have  also  calculated  the  losses  of  the  SRRs  as 
their  size  decreases.  It  is  found  that  losses  increase  as  the  size  of  the  SRRs  decreases.  We  have  also 
shown  that  the  effective  magnetic  permeability  ju{co)  does  not  approach  1  as  ft)  ^  0  and  co^oo.  This  is 
due  to  the  diamagnetic  response  of  the  SRRs.  We  present  numerical  results  which  support  this  conclu¬ 
sion.  One  way  to  increase  the  magnetic  resonance  frequency  co^  is  to  introduce  four-cut  SRRs  [20]. 
We  have  systematically  studied  co^  versus  the  size  of  the  four-cut  SRRs.  As  in  the  case  of  the  one-cut 
SRRs,  scaling  ~  ML)  is  obeyed  for  Z  >  500  nm  while  for  Z  <  500  nm  saturation  of  co^  is  also  ob¬ 
tained.  Using  the  retrieval  procedure  [21,  22]  we  have  calculated  ju{co),  both  the  real  and  imaginary  part 
of  ju{(d),  for  the  one-cut  and  four-cut  SRRs.  As  the  size  of  the  SRR  (for  both  the  one-cut  and  four-cut 
SRRs)  decreases  the  resonance  of  ju{co)  gets  weaker  and  the  imaginary  part  of  ju{co)  gets  stronger.  This 
also  means  that  the  losses  become  larger  as  the  size  of  the  SRR  gets  smaller.  The  paper  is  organized  as 
follows:  In  Section  II  we  present  the  results  for  one-cut  SRRs.  In  Section  III  we  present  the  results  for 
four-cut  SRRs.  Finally  in  Section  IV  we  give  our  conclusions. 


2  Scaling,  saturation  of  (o^  of  the  one-cut  SRRs 
and  their  magnetic  response 

We  simulate  single-ring  SRRs  similar  to  the  experimental  samples  previously  employed  in  the  200  THz 
frequency  range  [13]  for  near  infrared  and  optical  frequencies  using  the  commercial  finite-element  based 
PDF  solver  COMSOL  Multiphysics.  The  propagation  direction  of  the  incident  EM  field  is  along  the 
plane  of  the  SRR,  with  the  electric  field  parallel  to  the  sides  of  the  SRR  that  do  not  have  a  cut,  while  the 
magnetic  field  is  perpendicular  to  the  plane  of  the  SRR.  In  Fig.  1,  we  show  the  results  of  our  simulations 
for  the  resonance  frequency  co^  versus  the  inverse  size  of  the  unit  cell. 

For  large  enough  sizes  (Z  >  500  nm)  co^  scales  as  1/Z,  while  for  smaller  SRR  unit  cells  (Z  <  400  nm) 
co^  saturates  [20]. 


Unit  cel!  size  L 


Fig.  1  (online  colour  at:  www.pss-b.com)  Scaling  of  the 
magnetic  resonance  frequency  of  the  single-cut  SRR  with  the 
inverse  unit  size  is  shown.  The  parameters  of  single-cut  SRR 
are  given  by  length  /  =  0.8Z,  width  /  =  O.IZ,  thickness  of  the 
metal  t  =  0.05Z  and  gap-size  g  =  0.02Z .  The  linear  scaling 
(D^  oc  \IL  at  low  frequencies  (large  unit  cells)  breaks  down  at 
around  100  THz. 


2007  WILEY- VCH  Verlag  GmbH  &  Co.  KGaA,  Weinheim 


www.pss-b.com 


Original 

Paper 


phys.  stat.  sol.  (b)  244,  No.  4  (2007) 


1183 


Wavel  ength  [p  m]  Wave!  ength  [\i  m  1 


Fig.  2  (online  colour  at:  www.pss-b.com)  Real  part  (a)  and  imaginary  part  (b)  of  the  resonant  effective 
permeability  ju{cd)  are  shown  for  a  single-gap  SRR  as  function  of  frequency  for  different  sizes  L  of  the 
SRR  unit  cell:  Z  =  100pm  (red),  10  pm  (green),  1  pm  (blue),  167  nm  (pink),  100  nm  (turquoise)  and 
50  nm  (black). 

The  saturation  of  co^  at  high  frequencies  is  due  to  the  extra  inductance,  L^,  that  the  electrons  acquire 
as  the  SRR  unit  cell  gets  smaller  [20].  scales  like  ~1/L,  while  the  geometric  inductance  due  to  the 
shape  of  the  SRR,  scales  -L  and  the  capacitance,  C,  contributed  by  the  cut  in  the  SRR  rings  behaves 
~1/Z.  So  the  magnetic  resonance  frequency  co^  =  1/^ {L^+L^)C  =  \l\j q  +  €2!^  .  Thus  as  Z^O,  co^ 
saturates  to  a  constant  value. 

In  Fig.  2,  we  present  the  real  and  the  imaginary  part  of  the  effective  magnetic  permeability  ju{co)  for 
the  different  sizes  of  the  unit  cell  of  single-cut  SRRs.  Notice  that  the  resonant  response  in  ju^co)  gets 
weaker  as  the  size  of  the  system  decreases.  Notice  also  that  the  resonant  response  of  the  single-cut  SRR 
approaches  unity  for  high  frequencies.  The  smallest  size  we  have  studied  numerically  is  Z  =  50  nm.  Also 
the  peak  in  the  imaginary  part  of  ju{cd)  gets  smaller,  and  its  width  increases  as  the  size  of  the  single-cut 
SRR  gets  smaller.  The  same  is  true  for  the  magnetic  loss  tangent,  which  is  given  as  Im  ju{co)IRq  ju{co) 
and  shown  in  Fig.  3. 


Wavelength  [fim] 


Fig.  3  (online  colour  at:  www.pss-b.com)  Mag¬ 
netic  loss  tangent  Im  (ju(co))/Rq  (ju(co))  is  shown 
for  a  single-gap  SRR  as  a  function  of  frequency 
for  different  sizes  Z  of  the  SRR  unit  cell: 
Z  =  100  pm  (red),  10  pm  (green),  1  pm  (blue), 
167  nm  (pink),  100  nm  (turquoise)  and  50  nm 
(black). 


Wavelength  [/iin] 


Fig.  4  (online  colour  at:  www.pss-b.com)  Loss 
power  density  per  incident  power  flow  density 
(Poynting  vector)  is  shown  across  the  magnetic 
resonance  for  the  different  single-cut  SRR  sizes 
form  Fig.  1. 


www.pss-b.com 
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Wavelength  [fini] 

1000  100  10 

Fig.  5  (online  colour  at:  www.pss-b.com)  Real  part  of  the 
effective  magnetic  permeability  Re  (ju(cd))  is  shown  around 
the  magnetic  resonance  for  a  fixed  SRR  size  of  10  pm  and 
different  widths  of  the  metallic  ring.  Note  the  diamagnetic 
response  of  the  metal  area  (left,  below  resonance)  and  the 
total  enclosed  ring  area  (right,  above  resonance). 


I  10 

Frequency  [THz] 

A  better  way  to  assess  the  role  of  the  losses  in  the  SRR  is  to  calculate  directly  the  loss  power  density, 
which  describes  the  ohmic  loss  in  the  metal  and  is  given  by  P  =  crE^  -h  (dD/dt)  E.ln  Fig.  4,  we  present 
the  loss  power  density  normalized  to  the  incident  Poynting  vector  S.  Notice  that  the  loss  density  is  get¬ 
ting  larger  as  the  size  of  the  SRR  decreases. 

Another  curious  effect  shown  in  Fig.  2  is  that  ju(co)  does  not  approach  unity  for  low  frequencies  or  for 
high  frequencies.  In  Fig.  5,  we  show  the  real  part  of  the  effective  permeability  ju(co)  for  three  different 
SRRs.  The  size,  Z  =  10  pm,  is  the  same  for  all  three  cases.  The  only  difference  is  the  width  of  the  SRR 
ring  metal.  At  low  frequencies  ju(co)  deviates  from  unity  and  as  the  width  of  the  SRR  metal  increases 
ju((jL>)  approaches  0.5.  However,  in  the  high  frequency  regions,  just  below  the  peak  at  20  THz  ju(co)  is 
independent  of  the  width  of  the  SRR  metal.  To  understand  this  behavior,  we  have  systematically  studied 
the  current  densities,  as  well  as  the  component  of  the  magnetic  field  perpendicular  to  the  SRR  plane  for 
the  single-cut  SRR  at  low  and  high  frequencies.  In  Fig.  6a,  we  plot  the  current  density  of  the  single-cut 
SRR  of  10  pm  size  at  300  GHz  (i.e.  below  the  resonance). 

Notice  that  the  current  density  is  non-zero  only  at  the  outer  and  inner  edge  of  the  circular  SRR.  The 
current  direction  at  the  outer  side  is  clockwise,  while  at  the  inner  side  of  the  SRR  it  is  counterclockwise. 
In  the  gap  region  of  the  SRR  the  current  density  is  in  opposite  directions.  There  is  no  current  circulating 
inside  the  metallic  ring  at  300  GHz.  This  is  the  diamagnetic  response  of  the  single-cut  SRR.  The  dia¬ 
magnetic  response  can  be  easily  seen  in  Fig.  7a,  where  we  plot  the  z -component  of  the  magnetic  field. 


Fig.  6  (online  colour  at:  www.pss-b.com)  Current  distribu¬ 
tion  inside  the  SRR  is  shown  for  the  diamagnetic  response 
below  the  magnetic  response  at  (a)  300  GHz  (SRR  size 
Z  =  10  pm);  and  (b)  30  THz  (SRR  size  Z  =  100  nm);  as  well 
as  for  the  resonant  response  of  (c)  3  THz  (Z  =  10  pm),  and 
(d)  116THz(Z=  100  nm). 
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li)  (b)  (cj 


Fig.  7  (online  colour  at:  www.pss-b.com)  Magnetic  field  component  perpendicular  to  the  SRR  plane  is 
shown  in  the  SRR  plane  for  a  frequency  of  (a)  300  GHz,  below  the  resonance;  (b)  3  THz,  at  the  reso¬ 
nance;  and  (c)  13  THz,  above  the  resonance.  Note  the  negative  magnetic  response  in  (b)  and  the  two  dif¬ 
ferent  limits  of  diamagnetic  response  in  (a)  and  (c).  The  size  of  the  SRR  is  L  =  10  pm.  Yellow  color  repre¬ 
sents  positive,  blue  negative  and  green  zero  intensity  field. 

.  Notice  that  (green  color)  inside  the  metallic  ring,  as  expected  from  the  exclusion  of  the  mag¬ 

netic  field  from  the  interior  of  the  metal  due  to  the  diamagnetic  response. 

In  Fig.  7c,  we  plot  B^  for  the  same  size  SRR  (Y  =  10  pm)  but  at  higher  frequency  (13  THz,  above  the 
resonance).  At  13  THz  ju{co)  -0.5,  independent  of  the  width  of  the  metallic  ring  of  the  SRR.  As  can  be 
seen  from  Fig.  7c,  =  0  not  only  at  the  metallic  ring  of  the  SRR  but  also  throughout  the  enclosed  area 

of  the  SRR.  So  the  magnetic  field  is  excluded  from  the  area  of  the  SRR,  and  this  is  the  reason  that 
ju{co)  -  0.5  independent  of  the  width  of  the  metallic  ring.  Figure  7b  shows  the  strong  negative  magnetic 
response  of  the  SRR  slightly  above  the  resonance  frequency  co^  (3  THz). 

In  Fig.  6b  we  show  the  current  distribution  responsible  for  the  diamagnetic  response  for  a  smaller  SRR 
with  Z  =  100  nm  at  30  THz  (below  the  resonance).  Here  the  current  extends  deep  into  the  metal  of  the 
ring.  Figure  6c  and  6d  show  the  current  distribution  directly  at  the  resonance  frequency  for  the  10  pm 
SRR  (at  3THz)  and  the  10  pm  SRR  (at  1 16THz),  respectively. 

For  the  low  frequency  case  the  resonant  current  (which  overshadows  the  much  weaker  diamagnetic 
currents)  only  flows  in  a  thin  layer  (skin-depth)  at  the  inner  edge  of  the  SRR. 

At  high  frequency  (or  small  SRR  size),  the  skin-depth  as  well  as  the  resonance  frequency  starts  to 
saturate  (due  to  the  finite  electron  mass),  which  leads  to  an  increasing  ratio  of  skin-depth  over  SRR  size, 
hence  a  deeper  penetration  of  the  current  into  the  metal.  At  116  THz  (Fig.  6d)  the  resonant  current  is 
already  homogeneously  distributed  over  the  whole  cross-section  of  the  SRR  ring. 


3  Four-cut  SRR  give  negative  //  at  optical  frequencies 

As  we  discussed  in  Ref  [20]  the  4-cut  SRR  retains  the  negative  Re  (//(&>))  region  for  higher  frequen¬ 
cies  (up  to  about  550  THz).  As  one  can  see  from  Fig.  8,  the  resonance  frequency,  of  the  4-cut  SRR 


Fig.  8  (online  colour  at:  www.pss-b.com)  Scaling  of  the 
magnetic  resonance  frequency  of  the  four-cut  SRR  with  the 
inverse  unit  size  is  shown.  The  parameters  of  four-cut  SRR 
are  given  by  length  /  =  0.914Z,  width  and  thickness  of  the 
metal  w  =  t  =  0.251  L  and  gap-size  g  =  0.05L.  The  linear 
scaling  oc  HL  at  low  frequencies  (large  unit  cells)  breaks 
down  at  around  300  THz  which  is  higher  than  that  shown  in 
Fig.  1.  The  blue  solid  circular  symbol  (right  side  of  the 
vertical  dotted  line)  indicate  no  //  <  0  reached  anymore. 
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Frequency  {THz) 


Fig.  9  (online  colour  at:  www.pss-b.com)  Real  part  (a)  and 
imaginary  part  (b)  of  the  resonant  effective  permeability  ju{(d) 
and  loss  tangent  (c)  Im  (//(^u))/Re  {ju{(d))  are  shown  for  a  four- 
gap  SRR  as  function  of  frequency  for  different  sizes  L  of  the 
SRR  unit  cell:  L-1^  nm  (black),  56  nm  (red),  49  nm  (green), 
35  nm  (blue),  28  nm  (orange)  and  21  nm  (magenta).  In  part  (a), 
the  magnitude  of  magnetic  resonance  reduces  as  L  scales  down, 
and  only  L  =  10  nm  (black  curve)  shows  a  negative  //  at  around 
520  THz.  All  the  frequencies  below  520  THz  show  negative  ju . 


moves  to  higher  frequencies  for  smaller  SRR  size  but  slower  as  would  be  expected  from  linear  scaling. 
Simultaneously,  the  strength  of  the  resonance,  and  hence,  the  amplitude  of  the  resonant  ju(co),  decreases. 
Notice  that  even  at  SRR  sizes  of  21  nm,  there  is  a  small  resonance  of  //,  which  is  very  weak  so  that  jU 
does  not  reach  negative  values  anymore.  However,  jU  is  still  showing  some  structure  and  is  not  equal  to 
one  independent  of  frequency. 

In  Fig.  9b,  the  imaginary  part  of  ju(cd)  is  shown.  Notice  that  as  the  size  of  the  unit  cell  decreases,  the 
width  of  the  half-maximum  of  Im  (ju(co))  increases,  which  means  that  the  losses  are  getting  larger  as  the 
size  of  the  4-cut  SRR  decreases.  In  addition  the  4-cut  design  is  favorable  for  more-dimensional  materials 
not  only  for  its  highest  attainable  magnetic  resonance  frequency  but  also  for  its  inherent  symmetry.  For  a 
given  structural  size,  the  resonance  frequency  co^  increases  with  the  number  of  cuts,  which  brings  the 
resonant  refractive  index  closer  in  frequency  to  the  edge  of  the  Brillouin  zone,  =  ncIcoL,  promot¬ 
ing  periodicity  artifacts  [20,  21].  With  the  4-cut  design  the  resonance  wavelength  is  comparable  to  the 
size  of  the  unit  cell  and  the  effective  medium  theory  (2  »  Z)  is  not  obeyed  anymore. 

4  Conclusions 

Using  detailed  numerical  simulations  we  have  systematically  studied  the  scaling  of  the  resonance  fre¬ 
quency  CD^  of  single-cut  and  four-cut  SRRs.  For  both  cases,  it  is  found  that  co^  scales  inversely  propor- 
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tional  to  the  size  L  of  the  unit  cell  of  the  SRR.  For  the  single-cut  SRR  the  linear  scaling  is  obeyed  up  to 
100  THz,  while  for  the  4-cut  SRR  the  scaling  is  obeyed  up  to  500  THz.  Above  these  frequencies,  cd^ 
saturates.  The  reason  for  the  saturation  is  due  to  the  increased  kinetic  inductance  of  the  electrons  (com¬ 
pared  to  the  geometric  inductance  of  the  SRR  ring).  We  found  negative  effective  permeability  up  to 
500  THz  in  the  4-cut  SRR  design.  We  have  shown  that  the  losses  increase  as  the  size  of  the  SRR  unit  cell 
decreases.  The  saturation  of  the  resonance  magnetic  frequency,  as  well  as  the  losses,  are  responsible  for 
the  disappearance  of  the  negative  magnetic  SRR  response.  Finally,  we  numerically  calculated  the  current 
density,  and  the  z-component  of  the  magnetic  field  for  SRRs  at  low  and  high  frequencies.  For  high  fre¬ 
quencies,  the  current  density  is  non-zero  in  the  whole  metallic  ring  of  the  SRRs.  This  diamagnetic  re¬ 
sponse  can  explain  the  derivation  of  ju{co)  from  unity  which  is  observed  at  low  and  high  frequencies. 
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